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ABSTRACT 


Estimating wind speed characteristics plays an essential role in designing a 
wind power plant at a selected location. In this study, the Weibull, gamma, 
and exponential distribution models were proposed to estimate and analyze 
the wind speed parameters and distribution functions. Real measured data 
were collected from Medan City, Indonesia. The scale and shape factors of 


wind distribution for three years data were calculated. The observed 
cumulative probability of the three models was compared to predicted wind 
Keywords: characteristics. The probability density function (PDF) and the cumulative 
density function (CDF) of wind speed were also analyzed. The results 
showed that the Weibull model was the best model to determine PDF, while 
the exponential model was the best model to determine CDF for the Medan 
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1. INTRODUCTION 

Wind energy is one of the cheapest renewable energy, and it is a very large source of energy around 
the world. The global potential of wind energy is about 10 TW, and it can provide the world electricity needs 
as much as 35% [1]. According to IRENA reports [2], the world total installed capacity of wind energy is 
622.70 GW in 2019. While Asia has reached 40% of the global installed capacity [3], Indonesia has a total 
installed capacity of 76 MW in 2019, and this value is only 0.03% of the total installed capacity in Asia. 

Estimating wind speed characteristics is very important to assess the potential of wind energy [4], 
evaluate the performance of wind energy systems [5], and design a power plant at a selected location. 
Moreover, wind energy potential varies directly to the wind speed, and the small changes in wind speed will 
affect the amount of energy produced by the wind power plant. For these purposes, an appropriate wind 
model is required to accurately analyze the wind data measured at the site [6] because the wind energy is 
much related to the site and geographic conditions [7]. Another method to extrapolate wind data precisely 1s 
by using a site-specific wind power law [8]. 

The distribution function is one of the most widely used techniques, and it 1s vital for the analysis of 
performance losses in complex systems in engineering. In general, the distribution function can be used to 
predict and represent a variety of different phenomena. Moreover, it can also be applied to analyze the 
probability of wind speed. Characterizing the wind speed distribution function becomes simple when the 
measured histogram is correctly fitted with the analytical data [9]. The wind distribution function at certain 
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points can be interpreted as relative likelihood values. The wind distribution function also provides the 
frequency of wind speed [10], [11]. Wind distribution functions, such as the probability density function 
(PDF) and the cumulative density function (CDF) of the wind speed, are very crucial to determine the 
technical feasibility and the economic viability of a power plant project [12]. 

The popular distribution function to estimate wind energy is the Weibull distribution function. 
Weibull distribution function consists of shape and scale parameters, which can be calculated using several 
models available in the literature [13]-[15]. The shape factor identifies the width of the wind speed 
distribution and determines the peak wind distribution in any region [16]. The scale factor identifies the 
abscissa scale of wind distribution and the condition of most of the wind potentials in a particular location 
[17]. Mohammed et al. [18] applied the Weibull distribution function to analyze statistics on the maximum 
monthly wind speed and the minimum monthly wind speed in Zagora, Morocco. Mahmood et al. [19] also 
analyzed wind characteristics for the Al-Salman site in Iraq using Weibull distribution. 

In the wind energy application, a two-parameter Weibull is widely used and accepted to estimate 
wind energy potential. The two-parameter Weibull is represented in the form of a computable and flexible 
mathematical function [20]. Modeling wind speed with a two-parameter Weibull yielded significant values as 
reported in [21]-[24]. Teimourian et al. [25] applied a two-parameter Weibull to investigate other feasible 
wind sites in Iran. Sumair et al. [26] used a two-parameter Weibull with energy pattern factor method to 
estimate Weibull parameters.and assess wind potential in South Punjab. Sumair et al. [27] also compared 
two-parameter Weibull, Rayleigh, and lognormal distributions to model the wind potential of the coastal belt 
of Pakistan. The two-parameter Weibull performed the most accurately, and the Rayleigh distribution was 
found to be the least. 

Mert and C. Karakus [28] used the two-parameter Weibull, mixed Weibull, generalized gamma, and 
Burr 4-parameter distribution to describe the wind speed profile. They reported that the mixed Weibull 
distribution had a better match compared to the two-parameter Weibull [29]. In addition, estimating wind 
performance and assessment of the wind energy potential in wind energy planning projects using the two- 
parameter Weibull was not suitable for large-scale wind energy generation [10]. Three-parameter generalized 
gamma distribution [30] and normal distribution [31] were found to more accurately describe the wind 
characteristics compared to the two-parameter Weibull at different locations. Kollu et al. [32] developed a 
wind speed model based on Rayleigh and Weibull distributions. However, Weibull empirical data showed 
low confidence in the results, which prompted researchers to suggest various alternative analytic 
distributions, such as Gaussian inverse [33], lognormal [34], normal square root [35], Weibull inverse [36], 
the principle of maximum entropy [37], and Weibull-gamma composite [38]. Arreyndip et al. [39] suggested 
the general extreme distribution to analyzed wind energy variation and potential. In recent years several 
studies have been conducted to assess the potential for wind energy potential by using single and mixed 
distributions [35], [36]. 

It is obvious that no particular model fits all cases. For different wind locations, investigations 
should be carried out to obtain a suitable model for designing a wind power plant [40]. The Weibull 
distribution was the most used to determine the PDF of wind speed. The gamma and exponential distribution 
were often employed to determine the PDF of wind speed. However, there was still not enough information 
available about which model among the Weibull, gamma, and exponential distributions, is suitable to 
estimate PDF and CDF of the wind speed. This study aimed to propose and compare the three 
aforementioned models in estimating the best distribution functions of the wind speed characteristics in 
Medan City, Indonesia. 


2. RESEARCH METHOD 

In this study, the wind speed data was collected at Medan City (2°27°-2°47° М, 98°35’-98°44’ E), 
Indonesia. The historical monthly wind speed data for 2017-2019 was used to estimate the wind distribution 
function. Figure 1 shows daily wind speed from January 2017 to December 2019. Medan City is located at 
2.5-37.5 meters above sea level and has significant wind potential. Medan City has a tropical rainforest 
climate with dry and wet seasons that are not very distinct. Medan City has the driest months from February 
to April and experiences the wettest month in October [41]. The average wind speed in Medan City is around 
3.16 m/s with a standard deviation of 0.30433 [42]. 


2.1. Wind distribution functions 

Three distribution function models were proposed to evaluate wind speed in Medan City, namely 
the Weibull, gamma, and exponential distribution functions. Each distribution function has its own 
probability distribution function (PDF) and cumulative distribution function (CDF). The PDF is used to 
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determine the probability of a random variable falling within a specific range of values. In contrast, the CDF 
is used to determine the distribution of multivariate random variables. 
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Figure 1. Daily wind speed at Medan City for three years 


Weibull distribution function has the PDF as follows: 
k еру K-21 vik 
fv) ==() exp |-*] ,>k>0,v>0,c>0 (1) 
The corresponding CDF of the Weibull distribution function is defined by 
vik 
F(v) 21-exp|-?| , k>0,v>0,c>0 (2) 


where c and k are the scale and shape factors, respectively, and v is the observed wind speed (m/s). 

Like the Weibull distribution, the gamma distribution function also has two parameters, 1.е., a scale 
factor с > 0 and a shape factor k > 0 to control the wind speed distribution. The gamma PDF and CDF at 
observed wind speed v are defined by: 


f (v) = zz v" exp |-?]. k>0,v>0,c>0 (3) 
F(v) = f; 0. exp I-: dv, v >0,а є (0, оо) (4) 


where I (a) is the gamma function: 
l'(k) = f, v^ e-*dv, v0 (5) 
Different from the Weibull and gamma distributions, the exponential distribution function has only 


one parameter to control the wind speed distribution, 1.е., a scale factor c > 0. The exponential PDF and CDF 
at observed wind speed v are described as follows: 


f(v) = texp |-*], v>0,c>0 (6) 
F(v) =1-exp|-*|, v>0 (7) 


In a statistical method, the maximum likelihood (ML) can be applied to estimate the parameters of a 
probability distribution by maximizing a likelihood function [43]. According to this method, the shape (k) 
and scale (c) parameters are calculated to the total data points (№), as follows [44], [45]: 
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3. RESULTS AND DISCUSSION 

In this section, three proposed models, i.e., Weibull, gamma, and exponential, are applied to analyze 
wind speed parameters (scale and shape). The estimated wind speed distribution parameters of each model 
for the three years (2017-2019) are presented in Table 1, Table 2, and Table 3. Additionally, the yearly wind 
speed distribution parameters for 2017, 2018, and 2019 are also shown in the last row of Table 1, Table 2, 
and Table 3, respectively. 


Table 1. Estimated distribution parameters of wind speed for 2017 


Weibull distribution Gamma distribution Exponential distribution 
Scale (c) Shape (k) Scale (c) Shape (K) Scale (c) 
January 4.849 4.897 3.794 16.871 0.225 
February 5.343 4.667 3.593 17.560 0.205 
March 5.549 4.930 3.766 19.169 0.196 
April 5.049 5.729 5.280 24.674 0.214 
May SS 5.820 4.873 25.817 0.189 
June 5.834 4.783 3.820 20.391 0.187 
July 5.848 4.373 3.157 16.798 0.188 
August 5.848 5.056 4.194 21.324 0.197 
September 5.536 4.306 3.131 14.993 0.209 
October 5.261 4.631 3.690 17.484 0.211 
November 5.449 4.572 3.472 17.280 0.201 
December 4.585 4.556 3.525 14.767 0.239 
Yearly 5.373 4.560 4.356 21.618 0.202 


Table 2. Estimated distribution parameters of wind speed for 2018 


Weibull distribution Gamma distribution Exponential distribution 
Scale (c) Shape (k) Scale (c) Shape (k) Scale (c) 
January 4.849 4.897 3.794 16.871 0.225 
February 5.318 4.744 3.687 17.952 0.205 
March 5.573 4.959 3.835 19.601 0.196 
April 5.029 5.122 5.296 24.638 0.215 
May 9:797 5.857 4.985 26.601 0.187 
June 5.857 4.747 9:12 20.201 0.187 
July 5.867 4.341 3.099 16.537 0.187 
August 5.456 5.187 4.457 22.365 0.199 
September 5.354 4.161 2.905 14.128 0.206 
October 5.097 4.834 4.146 19.356 0.214 
November 5.490 4.540 3.455 17.312 0.200 
December 4.611 4.459 3.348 14.084 0.238 
Yearly 5.347 4.916 3.556 17.448 0.204 


Table 3. Estimated distribution parameters of wind speed for 2019 


Weibull distribution Gamma distribution Exponential distribution 
Scale (c) Shape (k) Scale (c) Shape (k) Scale (c) 
January 5.017 4.887 3.929 18.076 0.217 
February 5.389 4.053 2.805 13.702 0.205 
March 5.053 3.833 2.488 11.347 0.219 
April 5.228 5.303 4.395 21.146 0.208 
May 5.303 4.813 3.608 17.540 0.206 
June 5.866 3.122 3.867 20.862 0.185 
July 5.095 6.952 7.447] 35.495 0.210 
August 5.740 5.253 4.289 22.629 0.190 
September 5.702 6.714 6.532 34.787 0.188 
October 5.833 5.123 4.662 25.170 0.185 
November 5.596 5.688 4.833 25.036 0.193 
December 6.113 6.048 5.093 28.906 0.176 
Yearly 5.496 5.266 3.902 19.752 0.198 
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The monthly wind parameters of the three proposed models are also compared (as seen in Figure 2, 
Figure 3, and Figure 4). The parameters for 2019 are very different from other years, and this is due to 
climate change which changes the wind speed. Changes in wind speed are not followed by changes in wind 
parameters properly. Changes in wind speed parameters can result in changes in the wind distribution 
functions. As shown in Figure 2, Figure 3, and Figure 4, changes in parameters occur in March, May, July, 
and the following months, wherein these months have a significant difference in temperature in 2019 
compared to the other years. 

Furthermore, comparisons of the performances of three distribution models for 2017, 2018, and 
2019 by using the estimated parameters are shown in Figure 5, Figure 6, and Figure 7. The difference 
between the three models is very significant, especially in the exponential distribution model. This is due to 
the exponential distribution model only takes into account the scale factor, while the Weibull and gamma 
distributions use both scale and shape factors. 
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Figure 3. Estimated parameters of the gamma distribution for three years 
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Figure 4. Estimated parameters of the exponential distribution for three years 
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Figure7. Comparison of three distribution models for 2019 


The observed and expected cumulative probability for the three models in 2017, 2018, and 2019 are 
presented in Figure 8, Figure 9, and Figure 10, respectively. The expected cumulative probability of wind 
characteristics is a straight line and it is used as a data benchmark. Further investigations found that the 
exponential model is not around a straight line. This occurs because of the fluctuations in wind speed data 
due to climate change. 

From scale and shape factors and distribution probability results, the PDF and CDF of wind speed 
characteristics for Medan City are analyzed. For wind speed characteristic in 2017 (as seen in Figure 11), it is 
found that Weibull distribution PDF and CDF are 0.365 and 58%, respectively, with the wind speed of 5.2 
m/s. Gamma distribution gives PDF of 0.009 and CDF of 0.02%, with the wind speed of 7.2 m/s, while 
exponential distribution PDF and CDF are 0.045 and 74%, respectively, with the wind speed of 7.2 m/s. 

For wind speed characteristic in 2018 (as seen in Figure 12), the Weibull distribution PDF and CDF 
are 0.348 and 57%, respectively, with the wind speed of 5.2 m/s, while the gamma distribution gives PDF of 
0.000007 and CDF of 4%, with the wind speed of 7.2 m/s. In this case, the exponential distribution PDF and 
CDF are 0.046 and 77.5%, respectively, with the wind speed of 7.2 m/s. 

For wind speed characteristic in 2019 (as seen in Figure 13), the PDF and CDF obtained by the 
Weibull distribution model are 0.345 and 51%, respectively, with a wind speed of 5.2 m/s. In the same case, 
the gamma distribution gives PDF of 0.09 and CDF of 3%, with the wind speed of 7.2 m/s, while the 
exponential distribution gives PDF and CDF of 0.048 and 73.0%, respectively, with the wind speed of 7.2 
m/s. 

In addition, the PDF and CDF of the Weibull, gamma, and exponential distribution models are 
tabulated as in Table 4. Table 4 shows that the Weibull distribution is the best model to determine PDF, 
followed by the exponential and gamma distributions. On the other hand, the exponential distribution is the 
best model to determine CDF, followed by the Weibull and gamma distributions. 


Analysis of wind speed characteristics using different distribution models in Medan City... (Suwarno) 


1108 О ISSN: 2088-8694 


Weibull P-P Plot of Year 2017 Gamma Р-Р Plot of Year. 2017 Exponential P-P Plot of Year 2017 


Expected Cum Prob 
Expected Cum Prob 








of 0.0 0.2 0.4 06 08 10 1 
0.0 0.2 04 0.6 0.8 1.0 0.0 0.2 04 0.6 0.8 10 
Observed Cum Prob 
Observed Cum Prob Observed Cum Prob 





Figure 8. Observed and expected cumulative probability for the three models in 2017 
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Figure 9. Observed and expected cumulative probability for the three models in 2018 
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Figure 10. Observed and expected cumulative probability for the three models in 2019 


Table 4. Comparison PDF and CDF for three models 


Weibull model Gamma model Exponential model 

PDF CDF (%) PDF CDF (%) PDF CDF (%) 
2017 0.365 58.0 0.009 0.02 0.045 74.0 
2018 0.348 57.0 0.000007 4.00 0.046 Thà 
2019 0.345 51.0 0.09 3.00 0.048 73.0 
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Figure 11. PDF and CDF: (a) Weibull, (b) gamma, (c) exponential distribution models for 2017 
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Figure 12. PDF and CDF: (a) Weibull, (b) gamma, (c) exponential distribution models for 2018 
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Figure 13. PDF and CDF: (a) Weibull, (b) gamma, (c) exponential distribution models for 2019 


4. CONCLUSION 

Three different models have been applied to estimate parameters and distribution functions of wind 
speed in Medan City. Comparisons of the performances of three distribution models for 2017, 2018, and 
2019 wind data showed that the exponential model performed very differently from the Weibull and gamma 
models. The observed cumulative probability of wind characteristics by using the exponential model did not 
fit the expected values. However, the exponential distribution was found as the best model to determine CDF, 
followed by the Weibull and gamma distributions. On the other hand, the Weibull distribution was the best 
model to determine PDF. The gamma model gave the worst performances in determining both CDF and PDF 
values. 
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